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Abstract 

A supersymmetric construction of potentials describing the hard core inter- 
action of the neutron-proton system for low energies is proposed. It considers 
only the binding energy case and uses the approximation of the Yukawa po- 
tential given by Hulthen. Recent experimental data for the binding energy of 
the deuteron are used to give the involved orders of magnitude. 
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The neutron-proton (n-p) system is the simplest of the composite nuclei. It can 
exist in a stable bound state, the deuteron. This is a nuclide for which no excited 
bound states are known, i.e., it exists only in its ground state. For a n-p system in 
which there are no electrostatic forces (i.e., the force due to the two magnetic mo- 
ments can be considered as small as an irrelevant correction to the nuclear forces) the 
interaction is depicted by a short range central force. Forms of the related potentials 
commonly used are the square well, gaussian and Yukawa potentials (see for exam- 
ple [0].) Among them, the gaussian and the square well potentials have deserved 
special attention because their mathematical simplicity. On the other hand, the 
Yukawa potential has a deeper theoretical significance because its close connection 
with the 'exchange forces' responsible of the binding between nucleons. Historically, 
this potential is often quoted by the prediction of the existence of the n meson. 
Hulthen noticed that the Yukawa potential can be conveniently approximated by 
the expression 



where r is the distance between the nucleons, a > is a fixed length (the range of 
the potential) and Vo > is a fixed energy (the strength of the potential). 

For low energies, the Schrodinger equation for the Yukawa potential is not solv- 
able analytically and the potential ([!]) is used as the stand in for it (see || and 
references quoted therein). In this range of energies the problem of describing the 
n-p interactions can be characterized in two forms: that concerning to bound states 
or to scattering states. In the following, we shall only focus on the former case. 
On the other hand, empirically is known that the deuteron has a total spin 1. The 
obvious interpretation is that the spins of the neutron and the proton are parallel 
in the ground state, which would thus be described by a 3 Si-state. Therefore, we 
will restrict ourselves to the study of the nuclear binding energy ^-states. 

Let us stress now that the nucleons cannot approach each other closer than a 
certain distance, otherwise, the nucleus would not show an almost constant nuclear 
density. Then we must resort to phenomenological arguments [Q. The usual ap- 
proach is to modify the nuclear potential at small distances to be consistent with 
the experimental data. A simpler approach considers the so-called hard core model 
which is characterized by a short range infinite repulsion inside the attractive nu- 
clear interaction |5J. In other words, a realistic potential must contain more than a 
radial term of the form (|1|): it must also include an infinitely high barrier term. 

The present paper investigates the supersymmetric nature of the hard core 
deuteron's binding energy by doing calculations on the Hulthen's potential. We 
shall show that the new potential so derived presents a repulsive barrier term and 
can be either isospectral (e.g. it shares the same spectrum) or almost isospectral (the 
same spectrum except the ground state) to the Hulthen's potential. As a particular 
case, the new potential is chosen to be a representation of the hard core nuclear 
force describing the n-p low energy binding interactions. The involved orders of 
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magnitude are given by using the experimental data of —2.22456614(41) MeV for 
the deuteron binding energy £ d recently reported in |J . 

As usual for potentials depending on r, the Schrodinger equation for the Hulthen's 
potential ([!]) reduces to an eigenvalue equation for a particle in a one dimensional 
effective potential Ve(x) = £(£ + l)/x 2 — V(x), where £ is the azimuthal quantum 
number and x = r/a is a dimensionless radial coordinate. We are looking for the S 
states of binding energy (i.e., negative energies and £ = 0), therefore the following 
equation holds 

~ d 2 V 



+ u _ k 2 
dx 2 e x — 1 



i;(x) = 0, (2) 



where ip(x) = xR(x), with R(x) the standard radial wavefunction, and E = —k 2 = 
(2/za 2 /h 2 )£, Vo = (2fia 2 /h 2 )Vo, with \i the reduced mass. 

The standard procedure of solution carries out the eigenf unctions [0 
i/; n (x) = C n e- kx {l-e- x ) 2 F 1 (2k + l+n,l-n, 2k + l;e~ x ), n = l,2,..., (3) 
with C n a normalization constant 

The corresponding eigenvalues are given by 

E n = -hi = - (^^- ) , V >n 2 , n=l,2,... (4) 

Remark that the problem involves two mutually dependent parameters, Vq and k. 
In practice, the eigenvalue of the energy may be given by experiment while the 
strength Vo of the potential is to be determined. Therefore, equation ((4[) can be 
used to evaluate Vo in terms of E n . 

The dotted curve on Figure 1 represents the potential ([!]) allowing only one 
bound state, the deuteron's ground state, labeled by a subscript H. In Figure 2 we 
have plotted the well known corresponding probability density \iPh{x)\ 2 - 

As regards the supersymmetric scheme, we have in the first place to look for a 
superpotential w(x) solving the Riccati equation 

- w'(x) + w 2 (x) = V(x) - e, (5) 

where the prime denotes derivative with respect to x and the factorization energy 
e = —k 2 is, in principle, any real number. By a simple calculation we get for the 
particular solution 

w(x) = k-— -, k > 0, (6) 
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which is useful in the cases when the strength can be rewritten as Vq = 1 + 2k. As 
usual, the susy partner V of V is given by the shape invariance condition || 

V(x) = V(x) +2w'(x) } (7) 

leading to 

V{X) = + 2sinh'(x/2)- (8) 

Potential (|]) is a well known result in susy quantum mechanics (see |§). It has been 
used to study susy phase-equivalent potentials [[UJ and to establish some interesting 



connections between the susy and the variational method [11 



Observe now the appearance of the r.h.s. term in @. This term presents a 
singularity of order 2/x 2 at origin and behaves just as a repulsive centrifugal term 
with I = 1 in the neighborhood of x = 0. 

V( x )~-l±^ + l, x<<1 (g) 
x x l 

As the value x — 1 implies r = a, the approximation (^) holds in the range of the 
initial potential V(x). In the region x > 1, the potential V(x) rapidly becomes 
negligible (see Figure 1.) Here, the strength Vq = 1 + 2k plays the role of a coupling 
constant. 

The above results can be used to determine the eigenf unctions and eigenvalues 
connected with the new potential V(x). The procedure consits now in factorizing 
the related Hamiltonias by 

H = A^A + e, H = AA^ + e, (10) 

with 

A=4- + w(x). (11) 
ax 

It is straightforward to check that equations @, (|?]) and (|TT|) automatically lead to 
(|T0"|). Then, an intertwining relationship holds: 

HA = AH. (12) 

If ip(x) is an eigenfunction of H with eigenvalue E, equation (|l"2"D gives 

H(Aip) = E(Aip), Aip ^ 0. 

Therefore, if ip £ L 2 (R), we get the normalized eigenstate of 

= (E-e)^ 1/2 A^(x). (13) 

Now, let us stress on the information displayed by equations (|l0l)-(|i"3|). First, in the 
case when e = E n , for any n = 1, 2, . . ., the l.h.s. equation in fllCf ) applies on ipn{x) 
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as Hif} n (x) = E n ip n (x), and consequently A^Ai]) n (x) = 0. It is easy to check that 
Aip n (x) = is a sufficient condition to get square integrable functions. Therefore, 
equation (|13D means that i/j n ( x ) has not a susy partner ip n (x) and the couple of 
Hamiltonians (|10"D corresponds to of unbroken supersymmetry M. 

On the other hand, when e 7^ E n for every n = 1,2,..., there is no square 
integrable eigenfunction of H annihilated by A, and equation ( jlBf ) means that every 
i/j n {x) will have a susy partner ip n (x). Now, from the r.h.s. of (|i~0D, it is clear that a 
function ip*(x), obeying Aty* = 0, leads to Hip^x) = eip*(x). Hence, if ip* G L 2 (R), 
it must be added to the new set {ip}. When if)*(x) is an square integrable function 
the Hamiltonians (|H]) present unbroken supersymmetry, otherwise they correspond 
to of broken supersymmetry |§ . 

For the superpotential we are dealing with, one gets ip*{ x ) ^ e KX (l— e~ x )~ 1 , which 
is obviously not square integrable in [0, 00) for k > 0. Then, the supersymmetric 
behaviour of the Hamiltonians ( |i~0"|) lies in the selection of e, i.e., if e is chosen to be 
either an eigenvalue of H or not. 

In the following we shall consider the case when the initial potential V(x) allows 
the binding of only two states. The purpose of this convention will be apparent 
in the sequel. To get a system with only two energy levels, the strength of the 
Hulthen's potential has to be in the domain 4 < Vo < 9. Providing these values of 
V we get for the energies: —16 < E\ < —2.25, and —1.56 < E 2 < 0. 

In order to get an idea of the orders of magnitude involved, we note that, although 
there is no a priori reason why a should not be different for different sorts of the 
stationary systems described by V(x), we can take its numerical value as a = 3f. 
Such assertion is justified by the fact that the mean distance between nucleons 
(i.e., the size of the nucleus) is in the range of 2 or 4 / [0. Therefore, we get 
(h 2 /a 2 m p ) ~ 4.6113 MeV. Here, we have assumed that the neutron and proton 
masses are equal to m p , hence 2/i = m p . In this way, the experimental value of 
the deuteron binding energy Sd becomes in a dimensionless value — —0.4825 
{k d ~ 0.6946). Remark that, for the above selected domain of V , the deuteron 
energy lies in the domain of the exited state E 2 and not in the domain of the ground 
state Ei. 

Let us consider now the factorization energy fixed as e = Ei = —k\. In this 
case, as discussed above, we will have a couple (p~Q|) with unbroken supersymmetry. 
Hence Ai\)\{x) = 0, and the potential V(x), with Vq = 1 + 2k — 1 + 2ki, misses the 
ground state of V(x) and admits only one bound state (see equation (|T3|)): 

${x) = (E 2 - EiY 1 ' 2 

with eigenvalue 



— Inipzix) + w(x) 
ax 



(14) 
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We go a steep further and impose that the numerical value of E2 be determined 
by experiment and let it be equal to E d , the dimensionless value of the deuteron 
binding energy, therefore 

k t = Akd + 3 ~ 2.8892, V = Ak d + 4 ~ 6.7784 (16) 

which agrees with the previously stated domains for Vq, E\ and E 2 . Then, an initial 
potential ([[J), with range a = 3f and strength Vo — 31.2572 MeV, has a susy partner 
flU) allowing the binding of a single state if)(x) with an energy exactly equal to £ d - 
The main characteristic of this new potential is its centrifugal term, which reduces 
the range of approaching between the nucleons. Therefore, we have constructed a 
radial potential describing the hard core binding interaction between the nucleons. 

The probability density connected with the ground state (and single!) eigen- 
function ( |T4D of the hard core Hamiltonian H, for the numerical data mentioned 
above, has been plotted on Figure 2. Its features can be contrasted with those of 
the no core Hamiltonian H, labeled by |^#(x)| 2 , in the same figure. Observe the 
displacement to the right of ip with respect to ipn- An easy calculation shows that 
i/j(x) ~ (1 — e~ x ) iPh{%), hence, near the origin, %p goes to zero as x 2 whereas ipH 
goes just as x, and the probability to find the state ip(x) near to zero is minor than 
the probability to find i^h(x) at the same place. 

We shall now discuss some of the various implications of our results. First, the 
susy procedure accomplishes the derivation of a short range potential V(x) allowing 
only one bound state. This potential exhibits many of the qualitative features 
concerned with a hard core potential. It behaves predominantly as an effective radial 
potential Vi(x), with £ — 1, near the origin (see equation (|S])) and goes negligible 
with the increasing of x in the region x > 1. 

Second, from Figure 2, it is clear that there is a considerable probability of finding 
the two nucleons at distances larger than a. Therefore, the nuclear force connected 
with V(x) plays a relevant role only in a weak way because the neutron and the 
proton are outside each other's range so much of the time. That is, of course, a 
well known feature in the behaviour of the deuteron's ground state fl], [12|]. It is 



interesting to remark on the smoothest of ip(x), this is a nodeless function, just as 
one might expect for the eigenfunction of a single stable bound state. Therefore, 
the wavefunction (|TJ) corresponds to the hard core ground state eigenfunction of 
the deuteron. 

In summary, the method developed in the paper relied on the assumption that 
it is possible to describe the interacting system of a proton plus a neutron by a 
Schrodinger equation. The method allowed the construction of a hard core poten- 
tial for the deuteron and the estimation of the potential energy necessary to give the 
observed deuteron's binding energy. It is remarkable that even very refined experi- 
ments at low energies do not suffice to determine more than the range and strength 
of the involved potential, leading the detailed shape completely indeterminate. 
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Observe that the hard core hypothesis makes the strength (|16|), of the potential 
([l]), increased by V ~ SV h, where V h = 2k d + 1. In general, the nuclear forces are 
quite complicated and the assumption of a pure 3 S'-state does not suffice to explain 
the deuteron quadrupole moment, this makes necessary to introduce a tensor force 
[|J. However, it is well known that the Hulthen's potential gives a good approxi- 
mation for the binding energy in the terms discussed at the very begining of the 
paper. On the other hand, as mentioned above, there are different ways to modify 
the nuclear potential in order to get a hard core model, the potential derived here 
could be tested by the nucleon-nucleon scattering approaches, where the hard core 
hypothesis seems to be compatible with the empirical data || . 

As regards the unbroken susy results, we have to say that recent susy treatments, 
involving general solutions to the Riccati equation d^J), have shown the way to in- 
troduce susy partner potentials sharing exactly the same spectra Rl3| . That is, we 
could obtain a new Hamiltonian H which, joined in a susy couple with H, present 
unbroken supersymmetry by extending the superpotential @ to a general solution 
of (El). Although that method has been sucssesfully applied in the study of some 



interesting potentials (see e.g. [|1J, [15|]), it is out of the present scope and will be 
given elsewhere. 

I would like express my gratitude to L.M. Nieto for his unending patience during 
the long discussions concerning this paper. I had learnt quite a lot of him and 
J. Negro during my scientific stay in Valladolid (Spain). I am also greatly indebted 
to Miss M. Lomeli for her constant help in typing the manuscript. 

This work has been partially supported by Junta de Castilla y Leon, Spain 
(project C02/97). The kind hospitality at Departamento de Fisica Teorica, Univer- 
sidad de Valladolid, Spain, is also aknowledged. 



References 

[1] M.A. Preston and R.K. Bhaduri, Structure of the nucleus, (Addison, Read- 
ing, Mass. 1975); J.M. Blatt and V.E. Weisskopf, Theoretical Nuclear Physics, 
(Springer, N.Y. 1979) 

[2] L. Hulthen, Phys. Rev. 61, 671 (1942); L. Rosenfeld, Nuclear Forces, (North- 
Holland, Amsterdam 1948) 

[3] C.S. Lam and Y.P Varshni, Phys. Rev. A 4, 1875 (1971) 

[4] L. Hulthen and M. Sugawara, The two-nucleon problem in Encyclopedia of 
Physics, edited by S. Fliigge (Springer, Berlin 1957) vol. 39, p. 1. 

[5] R. Jastrow, Phys. Rev. 81, 165 (1951). 

[6] E.G. Kessler et. al, Phys. Lett. A 255, 221 (1999) 



6 



[7] S. Fliigge, Practical Quantum Mechanics (Springer, Berlin 1971); O.L. De Lange 
and R.E. Raab, Operator Methods in Quantun Mechanics (Claredon Press, Ox- 
ford 1991) 

[8] F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251, 267 (1995) 

[9] U. Laha, C. Bhattacharyya, K. Roy and B. Talukdar, Phys. Rev. C 38, 558 
(1988) 

[10] B. Talukdar, U. Das, C. Bhattacharyya and PK. Bera, J. Phys. A 25, 4073 
(1992). 

[11] E. Drigo Filho and R.M. Ricotta, Mod. Phys. Lett. A 10, 1613 (1995) 

[12] J.M. Levi-Leblond and F. Balibar, Quantics (North-Holland, Amsterdam 1990) 

[13] B. Mielnik, J. Math. Phys. 25, 3387 (1984); D.J. Fernandez C, Lett. Math. 
Phys. 8, 337 (1984); M.M. Nieto, Phys. Lett. B 145, 208 (1984) 

[14] D.J. Fernandez, L.M. Glasser and L.M. Nieto, Phys. Letts A 240, 15 (1998); 
D.J. Fernandez C. and V. Hussin, J. Phys. A 32, 3603 (1999) 

[15] J.O. Rosas-Ortiz, J. Phys. A 31, L507 (1998); J. Phys. A 31, 10163 (1998). 



7 



3 



6 




Figure 1: The susy partner potentials ([[]) and (||) with a = 3f 
and Vo — 31.2572 MeV. The deuteron's binding energy Sd — 
— 2. 2245 Me V has been ticked on the frame. The potential Vh has 
a strength Vh — 11.0173 MeV and the same value of a. 
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Figure 2: The deuteron's ground state probability density for 
the cases without core \iPh{x)\ 2 , and hard core |^(^)| 2 , with if)(x) 
given in ( T^) and the values displayed in Figure 1. Observe the 
displacement to the right of (^(a;)! 2 with respect to \iPh{x)\ 2 - 
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